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We briefly review the current status of the hadronic light-by-light scattering correction to the muon g — 2. 
Then we present our semi-analytical evaluation of the pion-pole contribution, using a description of the pion- 
photon- photon transition form factor based on large- Nc and short-distance properties of QCD. We derive a 
two-dimensional integral representation which allows to separate the generic features from the model dependence, 
in order to better control the latter. Finally, we sketch an effective field theory approach to hadronic light-by-light 
scattering which yields the leading logarithmic terms that are enhanced by a factor Nc- It also shows that the 
modeling of hadronic light-by-light scattering by a constituent quark loop is not consistent with QCD. 



1. Introduction 

The muon g — 2 is an important quantity that 
provides a stringent test of the Standard Model 
and which is potentially sensitive to new physics. 
However, for this purpose one first needs to well 
understand the hadronic contributions, i.e. vac- 
uum polarization effects and light-by-light scat- 
tering. The present picture of hadronic light-by- 
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Figure 1. The hadronic light-by- light scattering 
contribution to the muon g — 2. 
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light scattering is shown in Fig. 1 and the corre- 
sponding contributions to a M are listed in Table 1, 
taking into account the corrections made in the 
two full evaluations [1,2], after we had discovered 
the sign error in the pion-pole contribution [3,4]. 



Table 1 

Contributions to a M (xl0 10 ) according to Fig. 1. 
The last column gives the result when no form 
factors are used in the couplings to the photons. 



Type 


Rcf. [1] 


Rcf. [2] 


Ref. [3] 




(b) 


-0.5(0.8) 


-1.9(1.3) 




-4.5 


(c) 


8.3(0.6) 


8.5(1.3) 


8.3(1.2) 


+oo 




0.174 a 


-0.4(0.3) 






(d) 


1.0(1.1) 


2.1(0.3) 




~ 6 


Total 


9.0(1.5) 


8.3(3.2) 


8(4) b 





a Only a\ exchange. 

b Our estimate, using Refs. [1-3]. 



There are three classes of contributions to the 
hadronic four-point function [Fig. 1(a)], which 
can be understood from an effective field the- 
ory (EFT) analysis of hadronic light-by- light scat- 
tering [5,4]: (1) a charged pion loop [Fig. 1(b)], 
where the coupling to photons is dressed by some 
form factor (p-meson exchange, e.g. via vector 
meson dominance (VMD)), (2) the pseudoscalar 
pole diagrams [Fig. 1(c)] together with the ex- 



2 



change of heavier resonances (/ , en, . . .) and, fi- 
nally, (3) the irreducible part of the four-point 
function which was modeled in Refs. [1,2] by a 
constituent quark loop dressed again with VMD 
form factors [Fig. 1(d)]. The latter can be viewed 
as a local contribution -ipa^ipF^ to a M . The two 
groups [1,2] used similar, but not identical models 
which explains the slightly different results for the 
dressed charged pion and the dressed constituent 
quark loop, although their sum seems to cancel 
to a large extent and the final result is essentially 
given by the pseudoscalar exchange diagrams. We 
take the difference of the results as indication of 
the error due to the model dependence. 

Our approach to this problem consists of mak- 
ing an ansatz for the relevant Green's functions 
in the framework of large- Nc QCD. In this limit, 
an infinite set of narrow resonance states con- 
tributes in each channel. Then we perform a 
matching of the ansatz with chiral perturbation 
theory (ChPT) at low energies and with the oper- 
ator product expansion (OPE) at high momenta 
in order to reduce the model dependence. In 
practice, it is sufficient to keep a few resonance 
states to reproduce the leading behavior in ChPT 
and the OPE. In this way we show in Section 2 
that the pseudoscalar contribution now seems un- 
der control, due to our semi-analytical calcula- 
tion [3], using a pion-photon-photon form factor 
J r 7r o 7 * 7 * which fulfills the relevant QCD short- 
distance constraints [6], in contrast to the form 
factors used in Refs. [1,2]. 2 These findings are 
also corroborated by an EFT and large- Nc anal- 
ysis [4] which allows to calculate the leading and 
next-to- leading logarithms in a LbyL;had (g ec _ 3) 

2. Pion-pole contribution 

The contribution from the neutral pion inter- 
mediate state is given by the following two-loop 
integral (see Ref. [3] for all the details) 



X 



' d 4 qi f d 4 q 2 



1 



(2tt)V (2tt)4 q 2 q 2 ( qi +q 2 ) 2 
1 



[(j> + <7i) 2 _ m 2 ][(p — q 2 ) 2 - m 2 



2 Furthermore, the calculations in Refs. [1,2] were based 
purely on numerical approaches. 
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that involves the convolution of two pion-photon- 
photon transition form factors, see Fig. 1(c). The 
Tj are polynomials of up to sixth order in the 
momenta p, q\ , and q 2 , with p 2 = m 2 . 

Since no data on the doubly off-shell form fac- 
tor J r 7r o 7 » 7 , (q 2 , q 2 ) is available, one has to resort 
to models. In order to proceed with the ana- 
lytical evaluation of the two- loop integrals, we 
considered a certain class of form factors which 
includes the ones based on large-iVc QCD that 
we studied in Ref. [6]. For comparison, we have 
also used a vector meson dominance (VMD) and 
a constant form factor, derived from the Wess- 
Zumino-Witten (WZW) term. 

In large-TVp QCD, the pion-photon-photon 
form factor is described by a sum over an infinite 
set of narrow vector resonances, involving arbi- 
trary couplings, although there are constraints 
at long and short distances. The normaliza- 
tion is given by the WZW term, ^ r 7r o 7 , 7 , (0, 0) = 
-N c /(12n 2 F n ), whereas the OPE tells us that 

,{\ 2 q\(p-\qf) 



lim ^ r 7r o 7 . 7 . 



2F, 
3 q 



2 \a : 



+ 



A 2 A 3 
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In the following, we consider the form factors that 
are obtained by truncation of the infinite sum in 
large- Nc QCD to one (lowest meson dominance, 
LMD), and two (LMD+V), vector resonances per 
channel, respectively: 



^s* 7 D . r (??.9i)=- 



q\ 
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cv 



3 (g 2 



M 2 )' 



(3) 



(g?-g 2 2 ) = f { (dhM+qD+hMl+qlf 



+h 2 qjq 2 + h h {ql+q 2 2 ) + h 7 ^j j ((qj - M 2 J 



x(ql - M* 2 )(q 2 2 - M^iq 2 . 



(4) 
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with the constants cy = N c M v /(4tt 2 F 2 ) and the moduli of the Euclidean momenta 

h 7 = -N c M Vi M v J{Ait 2 F 2 ). The parameters 



hi,]%2, and /15 in the LMD+V form factor are 
not fixed by the normalization and the leading 
term in the OPE. We have determined these co- 
efficients phenomenologically [6,3]. In particular, 
jF x o 7>7 » (— Q 2 ,0) with one photon on-shcll behaves 
like 1/Q 2 for large spacelike momenta, Q 2 = —q 2 . 
Whereas the LMD form factor does not have 
such a behavior, it can be reproduced with the 
LMD+V ansatz, provided that hi = 0. A fit to the 
data yields moreover /i 5 = 6.93 ± 0.26 GeV 4 . An- 
alyzing the experimental data for the decay n° — > 
e + e~ leads to the loose bound |/i2|^20 GeV 2 . 

Note that the usual VMD form factor 
^ t {qUD ~ V[(«Zi 2 - M 2 )(q 2 My)] does 
not correctly reproduce the OPE in Eq. (2). 

For the form factors discussed above one can 
perform all angular integrations in the two-loop 
integrals analytically [3] using the method of 
Gegenbauer polynomials [7] . The key observation 
is that all form factors can be written as follows 



LbyL;, 



T^o Tr .{q 2 ,q 2 ) = f(q 2 ) - ^ 



~9M V% (q 2 ) 



4 



m 2 ; 



(5) 



This allows to cancel all dependences on q\ ■ q2 
in the numerators in J r 7r o 7 , 7 * (q 2 ,{qi+q2) 2 ) in the 
loop integrals in Eq. (1). Then one writes the 
propagators as (for Euclidean momenta K, L) 



1 



M 



(K — L) 2 +M 2 



I III n=Q 



with Zf L = (K 2 + L 2 + M 2 - [(K 2 + L 2 + M 2 ) 2 - 
4K 2 L 2 ] 1/2 )/(2\K\\L\), and uses the orthogonality 
properties of the Gegenbauer polynomials 



/ 



cM(K)C n (Q 1 -K)C m (K-Q 2 ) 



= 2ir z 



1 



C n (Q\ ■ Q2), 



(7) 



where for instance Q\ ■ K denotes the cosine of 
the angle between the four-dimensional vectors 
Qi and K . After performing the angular integra- 
tions in this way, the pion-exchange contribution 
to a,n can be written as a two-dimensional integral 
representation, where the integration runs over 



dQ 1 / alQ 2 
Jo 



x J>i(Qi,Q 2 ) /i(0i,Q 2 ), (8) 

i 

with universal [for the above class of form factors] 
weight functions u>i (rational functions, square 
roots and logarithms) [3]. The dependence on 
the form factors resides in /j. In this way we 
could separate the generic features of the pion- 
pole contribution from the model dependence and 
thus better control the latter. This is not possible 
anymore in the final analytical result (as a series 
expansion) for a^ yU '° derived in Rcf. [8]. Note 
that the analytical result has not the same status 
here as for instance in QED. One has to keep in 
mind that there is an intrinsic uncertainty in the 
form factor of 10 — 30 %, furthermore the VMD 
form factor used in that reference has the wrong 
high-energy behavior. 

The weight functions u>i in the main contribu- 
tion are positive and peaked around momenta of 
the order of 0.5 GeV. There is, however, a tail 
in one of these functions, which produces for the 
constant WZW form factor a divergence of the 
form ln 2 A for some UV-cutoff A. Other weight 
functions have positive and negative contribu- 
tions in the low-energy region, which lead to a 
strong cancellation in the corresponding integrals. 

In Table 2 we present the numerical results for 
the different form factors. All form factors lead to 
very similar results (apart from WZW). Judging 
from the shape of the weight functions described 



Table 2 

Results for a^ byL;7T " for the different form factors. 
In the WZW model we used a cutoff of 1 GeV in 
the divergent contribution. 



Form factor 


a L byL +< x 1Q 10 


WZW 


12.2 


VMD 


5.6 


LMD 


7.3 


LMD+V (h 2 = GeV 2 ) 


5.8 



4 



above, it seems more important to correctly re- 
produce the slope of the form factor at the origin 
and the available data at intermediate energies. 
On the other hand, the asymptotic behavior at 
large Qi seems not very relevant. The results 
for the LMD+V form factor are rather stable un- 
der the variation of the parameters, except for 
\i2- If all other parameters are kept fixed, our 
result changes in the range \h 2 \ < 20 GeV 2 by 
±0.9 x 10~ 10 from the value for h 2 = 0. 

Thus, with the LMD+V form factor, we get 



LbyL ; 7T 



+5.8 (1.0) x 10 



-10 



(9) 



where the error includes the variation of the pa- 
rameters and the intrinsic model dependence. A 
similar short-distance analysis in the framework 
of large- Nc QCD and including quark mass cor- 
rections for the form factors for the rj and rj' was 
beyond the scope of Ref . [3] . We therefore used 
VMD form factors fitted to the available data for 
^ r 7r o 7 » 7 * (— Q 2 , 0) to obtain our final estimate 

LbyL;PS „ LbyL; ,r , LbyL; I |_ ^LbyL;^' I 



= +8.3 (1.2) x 10" 



(10) 



An error of 15 % for the pseudoscalar pole contri- 
bution seems reasonable, since we impose many 
theoretical constraints from long and short dis- 
tances on the form factors. Furthermore, we use 
experimental information whenever available. 

3. EFT approach to a Jf yL;had 

In Ref. [4] we discussed an EFT approach to 
hadronic light-by-light scattering based on an 
effective Lagrangian that describes the physics 
of the Standard Model well below 1 GeV. It 
includes photons, light leptons, and the pseu- 
doscalar mesons and obeys chiral symmetry and 
U(l) gauge invariance. 

The leading contribution to a^ byLihad , of order 
p 6 , is given by a finite loop of charged pions with 
point-like electromagnetic vertices, see Fig. 1(b). 
Since this contribution involves a loop of hadrons, 
it is subleading in the large- Nc expansion. 

At order p s and at leading order in Nc, we 
encounter the divergent pion-pole contribution, 
diagrams (a) and (b) of Fig. 2, involving two 



WZW vertices. The diagram (c) is actually fi- 
nite. The divergences of the triangular subgraphs 
in the diagrams (a) and (b) are removed by in- 
serting the counterterm \ from the Lagrangian 
£( 6 ) = (a 2 /47r 2 F ) x V^Ts^tt + • • •, see the 
one-loop diagrams (d) and (e). Finally, there is 
an overall divergence of the two-loop diagrams (a) 
and (b) that is removed by a local counterterm, 
diagram (f). Since the EFT involves such a local 
contribution, we will not be able to give a precise 
numerical prediction for <+ byLihad . 




r 




A 



Figure 2. The graphs contributing to a^ yU,,a at 
lowest order in the effective field theory. 



Nevertheless, it is interesting to consider the 
leading and next-to-leading logarithms that are in 
addition enhanced by a factor Nc and which can 
be calculated using the renormalization group [4] . 
The EFT and large- Nc analysis tells us that 



+N C 




fM w ± M K ± 



1 2 . < Mo . 
In h ci In h c 



in. 



+ ^xlog'sj+ ^V c xlog'sj|, (11) 

where f(M^±/m fJ ,,M K ±/m fJ ,) = —0.038 represents 
the charged pion and kaon-loop that is formally 
of order one in the chiral and Nc counting and /i 
denotes some hadronic scale, e.g. M p . The coeffi- 
cient C of the log-square term in the second line is 
universal and of order Nc, since F 7T = 0(\/Nc)- 
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Unfortunately, although the logarithm is size- 
able, in a^ yUw0 there occurs a cancellation be- 
tween the log-square and the log-term. If we fit 
our result for the VMD form factor for large M p 
to an expression as given in Eq. (11), we obtain 



ln2 M, 



ci In — - + c 



= &c 

= (^) 3 C [3.94 - 3.30+ 1.08] 

= [12.3- 10.3 + 3.4] x 10~ 10 

= 5.4 x 10~ 10 , 



(12) 



which is confirmed by the analytical result of 
Ref. [8] (setting for simplicity M^o — m^): 

<CvmS° = [12 - 8.0 + 1.7] x 10- 10 = 5.7 x lO" 10 . 
This cancellation is now also visible in the re- 
vised version of Ref. [9]. In that paper the re- 
maining parts of ci have been calculated: c\ — 
-2 x (mo)/3 + 0.237 



-0.93±g;||, with our con- 



ventions for x an d x{M 



p Jcxp 



1.75 



+ 1.25 
-1.00- 



Finally, the EFT analysis shows that the mod- 
eling of hadronic light-by-light scattering by a 
constituent quark loop is not consistent with 
QCD. The latter has a priori nothing to do 
with the full quark loop in QCD which is dual 
to the corresponding contribution in terms of 
hadronic degrees of freedom. Equation (11) 
tells us that at leading order in Nq any model 
of QCD has to show the behavior a ^ byUh,id ~ 
(a/ir) 3 Nc[N c ml/(48ir 2 F%)} In 2 A, with a univer- 
sal coefficient, if one sends the cutoff A to infin- 
ity. From the analytical result given in Ref. [10] 
for the quark loop, one obtains the behavior 

a L byL; c Q M _ ( a /^)3 A T c ( m 2/ M 2 ) + ^ for Mq ^ 

to m , if we interpret the constituent quark mass 
Mq as a hadronic cutoff. Even though one may 
argue that N c / 4:8ir 2 F% can be replaced by 1/Mq, 
the log-square term is not correctly reproduced 
with this model. Therefore, the constituent quark 
model (CQM) cannot serve as a reliable descrip- 
tion for the dominant contribution to a Lb y L ; had^ 
in particular, its sign. Moreover, we note that 
the pion-pole contribution is infrared finite in the 
chiral limit 3 , whereas the quark loop shows an 
infrared divergence In (Mq /m M ) for Mq — > [10]. 

3 This can be shown by studying the low momentum be- 
havior of the weight functions Wi corresponding to the 



4. Conclusions 

The pseudoscalar pole contribution a J ^' yIj ' ps 
seems to be under control at the 15 % level. More- 
over, the EFT and large-A<c analysis provides a 
systematic approach to a Lb y L ; had an( j yj e y s the 
leading and next-to-leading logarithmic terms, 
enhanced by a factor Nc- It also shows that the 
modeling of hadronic light-by-light scattering by 
a constituent quark loop is not consistent with 
QCD. Since model calculations for the dressed 
charged pion and the dressed constituent quark 
loop yield slightly different results, our present 
estimate reads (by adding the errors linearly) 

a LbyL;had = + g ^ x JQ-10 _ ^3) 
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two-loop diagrams 2(a)-(c) and the one-loop diagrams 
2(d)+(e) (given in Ref. [4]) for M^a -> 0, see also Ref. [9]. 



